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1. Motivation and Objective 

1) Modeling of the near wall turbulence. We aim to develop a second order clo- 
sure for near wall turbulence. As a first step of this project, we try to develop a 
k - € model for near wall turbulence. We require the resulting model to be able to 
handle both near wall turbulence and turbulent flows away from the wall, compu- 
tationally robust, and applicable for complex flow situations, flow with separation, 

for cxfitinplc. 

2) Modeling of the bypass transition. We aim to develop a bypass transition 
model which contains the effect of intermittency. Thus, the model can be used 
for both the transitional boundary layers and the turbulent boundary layers. We 
require the resulting model to give a good prediction of momentum; and heat transfer 
within the transitional boundary and a good prediction of the effect of freestream 
turbulence on transitional boundary layers. 

2. Work Accomplished 

In the past year, progress has been made in both topics mentioned above (i.e., 
modeling of the near wall turbulence and modeling of the bypass transition). In t e 
paragraphs below, these two topics will be reported separately. 

2.1 Modeling of Near Wall Turbulence 

Because of the wide range of scales involved in a turbulent flow, DNS (direct 
numerical simulation) is limited to flows of moderate Reynolds number and simple 
geometry. Turbulence modeling is the only viable approach for the calculation of 
turbulent flows of engineering interest. In turbulence modeling, the k-e model is the 
most widely used model in engineering calculations. The Standard k-e Model ’ 
was devised for high Reynolds number turbulent flows and is traditionally used 
in conjunction with a wall function when it is applied to wall bounded turbulent 
flows. However universal wall functions do not exist in complex flows and it is thus 
necessary to develop a form of k-e model equations which can be integrated down 

to the wall. , , 

Jones and Launder 3 were the first to propose a low Reynolds number k e model 

for near wall turbulence, which was then followed by a number of similar k - e 
models A critical evaluation of the pre-1985 models was made by Patel etal. . More 
recently proposed models can be found in Shih 5 and Lang and Shih • Three major 
deficiencies can be pointed out about existing k-e models. (Some of the models may 
have only one or two of the three deficiencies.) First, a near wall pseudo-dissipation 
rate was introduced to remove the singularity in the dissipation equation at the 
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wall. The definition of the near wall pseudo-dissipation rate was quite arbitrary. 
Second, the model constants were different from those of the Standard k — e Model, 
making the near wall models less capable of handling flows containing both high 
Reynolds number turbulence and near wall turbulence, which often occurs for real 
flow situations. Patel et al. 4 put as the first criterion the ability of the near wall 
models to predict turbulent free shear flows. Third, the variable y + is used in the 
damping function of the eddy viscosity formulae. Since the definition of y + 
involves u T , the friction velocity, any model containing y + cannot be used in flows 
with separation. 

Effort is made to propose a new k—e model for near turbulence which is free of the 
three deficiencies mentioned above. In this model, k 1 ! 2 is chosen as the turbulent 
velocity scale. The time scale is bounded from below by the Kolmogorov time scale. 
The dissipation equation is reformulated using this time scale and no singularity 
exists at the wall. Thus, it is no longer necessary to introduce the near wall pseudo- 
dissipation rate. The model constants used are the same as in the Standard k-e 
Model. Thus, the proposed model will be also suitable for flows away from the wall. 
An earlier version of the model, which contains y + in the damping function, was 
proposed and reported in Yang and Shih 7 . The model is now improved by using 
Ry = k l l 2 y/v instead of y + in the damping function. Hence, the present model can 
be used for flows with separation. 

In the present model, the eddy viscosity is given by 

vt = c^fukT ( 1 ) 

where the time scale T is written as 

T=*+(^) 1/2 - (2) 

The first part is the time scale conventionally used for high Reynolds number tur- 
bulent flows and the second part is the Kolmogorov time scale. Away from the wall, 
the first part is much large than the second part while near the wall the second part 
dominates, giving the Kolmogorov time scale as the turbulent time scale at the wall. 
The time scale given is bounded from below by the Kolmogorov time scale and is 
always positive. 

The damping function / M is given by 


U = [ 1 — exp(-aiR y - a^Ry - a 5 Rj)] 1/2 (3) 

where a\ = 1.5 x 10 -4 , a 3 = 5.0 x 10~ 7 , a 5 = 1.0 x 10 -1 °. The damping function 
is chosen such that the shear stress has the correct near wall asymptotic behavior. 
Away from the wall, approaches one as required. 

The modeled transport equations for k and e are 

k -I- Ujk.j = \(v + — )k,j).j- < UiUj > Ui,j - c, (4) 
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e + Uje.j = [{v + + {-C u < W > V* ~ C ^ e )! T + "VrUijkUw (5) 

Because T is always positive, the dissipation rate equation does not have a singu- 
larity at the wall. , A 

Turbulent channel flows and flat plate boundary layer flows at different Reynolds 

numbers were calculated using the proposed model. The mean velocity, turbu- 
lent kinetic energy, turbulent shear stress and turbulent dissipation rate for tur- 
bulent channel flow at Re r = 395 and turbulent flat plate boundary layer flow at 
Re e = 1410 are shown in Fig 1 and Fig 2, respectively. DNS data for these cases 
are shown for comparison. Also shown are the predictions using the Jones-Launder 
model and the k - e model proposed by Chien*. These two models are chosen 
because the Jones-Launder model is the first k - e model for near wall turbulence 
while Chien’s model is known to perform quite well for turbulent boundary ayer 
flows Overall, the proposed model is found to give a better prediction. Calcu a- 
tions were also made for turbulent flat plate boundary layers at larger Reynolds 
numbers, turbulent boundary layers with pressure gradient (favorable pressure gra- 
dient, adverse pressure gradient, and increasingly adverse pressure.) The results of 
these computation and the comparisons with the available experimental data can 
be found in Yang and Shih 9 . 


2.2 Modeling of Bypass Transition 

In a quiescent environment, transition is preceeded by the amplification of Tollmien- 
Schlichting waves. These waves eventually break down, giving rise to turbulent 
spots, which can be viewed as the onset of transition. In an environment with high 
freestream turbulence, say the flow passing over a turbine blade, turbulent spots 
are formed due to the transport of turbulence from the freestream to the boundary 
layer rather than the T-S wave amplification. This type of transition is called bypass 
transition. Accurate prediction of bypass transitional boundary layers is very im- 
portant for internal fluid mechanics because a significant proportion of the turbine 
blade is in the a transitional boundary layer region. Furthermore, the performance 
and the life of a turbine are directly related to the peak vales of the momentum and 
heat transfer both of which occur in the transitional boundary layer. 

Priddin 10 was the first to notice that the low Reynolds number two equation 
models have the potential to predict transitional flows under the influence of the 
freestream turbulence. This is probably due to the fact that the generation o 
turbulent spots in a boundary layer is a random process and the flow is almost 
fully developed turbulent within a turbulent spot. A detailed calculation procedure 
was given by Rodi and Scheuerer 11 , in which the Lam & Bremhorst low Reynolds 
number k - e model was used. More recently, a comparative study of the perfor- 
mance of existing low Reynolds number k-e models in predicting laminar-turbulent 

transition was made by Fujisawa 12 . ... 

While the low Reynolds number k-e models could mimic transition, the quan- 
titative predictions do not compare very well with the experimental data. This is 
due to the fact that all these low Reynolds number k-e models were originally 
proposed for fully developed turbulent flows and did not take into consideration the 
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distinct feature of a transitional boundary layer - intermittency. The intermittency 
of a transitional boundary layer is measured by the intermittency factor which can 
be viewed as the percentage time a transitional boundary layer is in the turbulent 
state due to the passing of a turbulent spot. 

We propose a model for the calculation of transitional boundary layers, which 
takes the effect of intermittency into consideration. The model is based on the k-e 
model for near wall turbulence we have stated above. The effect of intermittency 
is introduced through the following argument: since the percentage of time that a 
transitional boundary layer is turbulent is measured by the intermittency factor, 
a model for transitional boundary layers could be constructed from a model for 
turbulent boundary layers by multiplying all the terms due to turbulence mechanism 
by a weighting factor 7, which is linearly related to the intermittency factor. Thus, 
the governing equations for the fiat plate transitional boundary layers are, after 
using the boundary layer approximation and the eddy viscosity assumption, 


dU 

dx 


dV n 
+ — 0, 
dy 


( 6 ) 


TT dU „dU d \, , , 

c 'aJ + v % = aSr + w) 


dU 

dy J ’ 
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dx dy dy [ 
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dU, 2 


+ 7M^) 2 -7e, 
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+ 7 
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1 ,dU ... 


where x , y are the coordinates along and normal to the plate and U, V are the mean 
velocities in the x, y directions, respectively. 

In order to close the above equations, an expression for the weighting factor 7 
is needed. The weighting factor is assumed to be related to both the freestream 
turbulent level and the intermittency factor of the boundary layer. The intermit- 
tency factor is assumed to be determined by the local state of the boundary layer. 
We use H, the shape factor, to characterize the local state of the boundary layer 
since both the intermittency factor and the shape factor change monotonically from 
the laminar boundary layer to the turbulent boundary layer. Experimental results 
by Abu-Ghannam and Shaw 13 are used as a guide to construct this function. The 
weighting factor is also assumed to change in the y direction in such a way that 
outside the boundary layer, the weighting factor is one since the freestream is gov- 
erned by decaying turbulence. The final expression for the weighting factor is given 
in Yang and Shih 14 . 

The above system of parabolic equations need to be supplymented by boundary 
conditions at the wall and at the freestream and by initial conditions at the starting 
point of the calculation. At the wall 


U = V = k = 0, 
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e = ' (10) 

€w 7 n dy 1 

The boundary condition for e was obtained by applying the equation for turbulent 
kinetic energy down to the wall. 

At the edge of the boundary layer, the flow variables are given by the freestream 

values, i.e. / 1t \ 

U = U e , k = fc e , e = e e . ( u ) 

The fc e e e in the above are found from the transport equations for k and e, with 
the condition that the gradient of the flow variables in the y direction vanishes as 
the free stream is approached. Thus, 


TT — = (13) 

e dx T 

kt0 and £«o (the values of K and c« at the leading edge, for example) ate needed. 
k M is obtained from the experiment, and e.o is determined in such a way that the 

resulting k e (x ) profile agrees with the experiment. , 

One of the issues in the calculation of transitional boundary layers through the 
low Reynolds number k - e models is the prescription of the initial profiles for 
the turbulence kinetic energy and its dissipation rate, the later of which could no 
be found from the experiment directly. An expression for the initial profiles were 
given in Rodi and Scheuerer 11 . However, computations by Yang and Shih which 
tested the effect of the initial conditions on the transition prediction found, in 
agreement with the findings of Patankar and Schmit 16 , that the predicted onset of 
the transition is sensitive to the initial profiles. This sensitivity of the results to the 
initial conditions suggests that the only place where the initial conditions could be 
specified unambiguously is at the leading edge. At the leading edge, the turbulent 
kinetic energy and its dissipation rate take constant profiles, the values of which are 

determined by the law for the decaying turbulence. 

With the initial conditions given at the leading edge and the boundary conditions 
riven above, the solutions are marched downstream. Flat plate boundary layers with 
free stream turbulence levels of 3% (Case T3A) and 6% (Case T3B) respectively 
were calculated using the present model. These are the benchmark cases m an 
ongoing project coordinated by Savill 17 , testing the capability of turbulence models 
in predicting transitional flows. Fig. 3 shows the variation of skin friction coefficient 
Cf against Re x . Results from the experiment is shown for comparison. In addition, 
the prediction of the Launder-Sharma model is also shown in the figure because it 
was reported that among the lower Reynolds number k - e models, the Launder- 
Sharma model performs best for transitional boundary layers. It is clear that the 
present model gives a better prediction. Other features in the transitional boundary 
layers and the calculations of the transitional boundary layers with other levels of 
freestream turbulence can be found in Ref. 14. 
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3. Future Plans 

1) Modeling of the near wall turbulence. 

The proposed k-e model is only tested for simple parabolic flows so far. Because 
of the form of the model equations, the proposed model can be used in complex 
flow situations, flow with separation for example. The performance of the proposed 
model in those situations will be tested. 

We will work on the second order closure for near wall turbulence. In particular, 
we will be looking at the effect of the mean flow inhomogeneity on the pressure 
strain correlation. We are hoping to represent this effect rationally, so that the 
ad hoc damping functions currently being used in all the near wall second order 
closures can be avoided. 

2) Modeling of bypass transition. 

We will apply the proposed model to transitional boundary layers with pressure 
gradient and curvature. We will also extend the model to thermal boundary layers. 
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Fig lc: Shear stress profile for channel flow at Re r = 395. 
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Fig 2c: Shear stress profile for channel flow at Reg — 1410. 
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Fig 2d: Dissipation rate profile for channel flow at Reg — 











